Electron Dynamics in a Coupled Quantum Point Contact 
Structure with a Local Magnetic Moment * 

Vadim I. Puller and Lev G. Mourokh 

Department of Physics and Engineering Physics, 
Stevens Institute of Technology, Hoboken, New Jersey 07030, USA 

A. Shailos and J. P. Bird 

Department of Electrical Engineering and Center for Solid State Electronics Research, 
Arizona State University, Tempe, Arizona, 85287-5706, USA 

Abstract 

We develop a theoretical model for the description of electron dynamics in coupled quantum 
wires when the local magnetic moment is formed in one of the wires. We employ a single-particle 
Hamiltonian that takes account of the specific geometry of potentials defining the structure as well 
as electron scattering on the local magnetic moment. The equations for the wave functions in both 
wires are derived and the approach for their solution is discussed. We determine the transmission 
coefficient and conductance of the wire having the local magnetic moment and show that our 
description reproduces the experimentally observed features. 
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I. INTRODUCTION 



The low-temperature conductance of quantum point contacts (QPCs) is well known to 
be quantized in units of 2e^/h, a phenomenon that can be explained in terms of a simple 
transmission (Landauer) picture in which the influence of electron-electron interactions is 
neglected jl|. While this model is remarkably successful in accounting for the observation 
of conductance steps at integer units of 2e^/h, it is unable to explain the origin of the 
additional conductance plateau, observed near 0.7 ■ 2e^ /h in numerous experiments. (For 
an overview of this issue, see Ref. ^j.) While many different theoretical models have been 
proposed to account for the origins of the 0.7 feature, there is a wide consensus that it 
should be associated with some novel many-body effect. In particular, there is growing 
consensus that this feature is associated with the development of a net magnetic moment 
m the QPC 0,y,3. In our recent work j^, we have explored the use of coupled quantum 
wires as a means for providing electrical detection of the local-moment formation. The 
device structure that we have studied is shown in Fig. 1 and was formed in the two- 
dimensional electron gas of a GaAs/AlGaAs quantum well. The device was realized by 
means of electron-beam lithography, and lift-off of Ti-Au gates. These gates were formed on a 
Hall bar with eight ohmic contacts, positioned uniformly along its upper and lower edges. In 
suitable combinations, these contacts could be used to make four-probe measurements of the 
conductance of either wire, or of the quantum dot itself (as indicated in Fig. 1). Of particular 
interest here is the non-local measurement (right panel) that can be made by measuring 
the conductance through one (fixed) wire as the gate voltage {Vg) applied to the other 
(swept) wire is varied. The key result of our experiment is that as the swept wire pinches 
off, a resonant enhancement of the conductance of the fixed wire is observed. A convincing 
theoretical explanation for this resonant interaction was provided in a subsequent theoretical 
report by our group \}\. Based on a modified Anderson Hamiltonian, we have shown that the 
resonant interaction with the local magnetic moment formed in the swept wire leads to an 
additional positive contribution to the density of states of the fixed wire and, consequently, to 
an enhancement of its conductance. While this analysis provides a qualitative understanding 
of the resonant interaction between the quantum wires, the tunnel matrix elements involved 
in the Anderson Hamiltonian are generally unknown and have thus far been used as fitting 
parameters. The influence of the specific device geometry on these matrix elements has 
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thus far been neglected, even though geometry-related effects are known to be important 
for the description of scattering in one-dimensional structures 0, ^. To overcome these 
shortcomings, in the present paper we develop a more comprehensive theory for electron 
dynamics in the coupled-wire system and attempt to calculate the amplitude of the resonant 
inter-wire interaction from first principles. The basic idea of this approach will be to calculate 
the single-electron transmission properties in a device potential that is modified by the 
presence of an extra scattering term, arising from the presence of a local magnetic moment in 
one of the wires. The formulation of this idea is given in Section II where we derive equations 
describing the dynamics of electrons in the swept and fixed wires. The approach to treat 
these equations is given in Section III, where the expression for conductance is obtained in 
terms of transmission matrix elements. In the present paper we restrict ourselves to the 
analysis of electron dynamics in the swept wire with the examination of the fixed wire to be 
published elsewhere jlOf, and in Section IV we determine the transmission coefficient and 
conductance for the swept wire. In this, we obtain experimentally observed features such as 
additional 0.75 ■ le^jh plateau for the ferromagnetic coupling and 0.25 ■ 2e^/h plateau for 
the antiferromagnetic coupling 11, 3|- The conclusions are presented in Section V. 



II. ELECTRON MODES IN THE COUPLED QUANTUM WIRE STRUCTURE 

We start our description of electron dynamics in the coupled quantum wire structure 
from the following single-particle Hamiltonian 

Ho = K, + Ky + U{x) + W{y) + V{x, y) - J{x, y)& ■ I, (1) 

where and Ky are the kinetic energy operators for the electron localized in 2D plane, 
W{y) is the double- well potential describing the two quantum wires (Fig.2, center panel), 
V{x,y) is the potential of the tunnelling channel connecting the two wires (Fig.2, right 
panel), and U{x) describes the smooth bottleneck shape of the quantum wire channels. 
The last term simulates exchange coupling between the conductance electrons (Pauli ma- 
trices (?) and the local moment, S, which is assumed to be a spin-1/2 magnetic mo- 
ment with J{x, y) as the coordinate-dependent exchange coupling constant. The potentials 
?7(x), J{x,y), and V{x,y) — > as x — > ±oo. The potential V{x,y) is very sharp in com- 
parison with the variation of U{x) in the x-direction. J{x,y) has an x-dependence similar 
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to that of U{x), since the spatial characteristics of the local magnetic moment formed in the 
conducting channel are determined by the shape of this channel. 
We write the Schrodinger equation in the form 



Ho^{x,y) = Eip{x,y), 



(2) 



where the symbol " hat" in this and other equations is used for operators and wave functions 
in the four-dimensional spin space of the two spins. The basis vectors in this space are given 



by 
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XI = lie) \h) , X2 = lie) lis) , X3 = |Te) \ls) , and ^4 = lie) IT5) , (3) 

where |te) (lie)) and \^s) ills)) are spin- up (spin-down) states of the electron spin, a, and 
the local moment spin, S, respectively. 

The solution of the Schrodinger equation, Eq. (j2)), can be expanded in terms of the spin 
functions, Eq. (jS)), as 

4 

^{x,y) = J2xa^a{x,y). (4) 

It should be noted that the Hamiltonian, Eq. (^, is not diagonal in the spin space 
determined by the uncoupled representation, Eq. Q, due to the presence of the exchange 
term. This term can be diagonalized by means of a unitary transformation to the coupled 
representation with transformation operator 

/1 n n n \ 
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(5) 



V 

The wave function in the coupled representation is given by 

i^'ix^y) = X~^ip{x,y) 



(6) 



with ipaicn = 1,2,3) describing the triplet and with ip'^ describing the singlet spin states. 
It should be emphasized that for x — > ±00 the Hamiltonian is diagonal in the uncoupled 
representation due to the vanishingof potentials U{x) = 0, V{x,y) and J{x,y) = 0. 

Following the procedure of Ref. ja| we expand the full wave functions in terms of different 
propagating modes 

i^ix,y) =Y,0n{x)<l>n{y) (7) 



with the transverse structure of n-th mode given by the solutions of the equation 

[Ky + W{y)]^{y) = E^^„{y). (8) 
Correspondingly, the wave functions ipn{,x) obey the coupled equations 

[E - En - K^- Un{x)\ (pn{x) = f Km(a;) - Jnm{x)a ■ Sj (pm{x) (9) 

where 

Vnm{x) = J dyKiy)V{x,y)<i>Uy), (10) 
Jnmix) = J dy<l>l{y)V{x,y)<l>Uy), (H) 

and Un = U{x) + Vnn{x). 

In the following analysis we make a number of simplifications in Eq. First, we note 
that if the wires are well separated, the wave functions are strongly localized in one 

of the two wires, therefore we can distinguish the modes propagating in each of the wires. 
We assume that the shape of the confining potential W{y) is such that one of the wires is 
close to pinch off {swept wire), i.e. it has only one propagating mode (described by the wave 
function 0o{x)) with the transverse confinement (subband bottom) energy. En, less than 
the Fermi energy, whereas the other wire {fixed wire) has several propagating modes. The 
localized magnetic moment is supposed to form in the only subband of the swept wire, hence 
the exchange coupling can be approximated as Jnm{x) = Sn,oSm,oJ{x)- Thus the system of 
equations is reduced to 



E-Eo-K^- Uo{x) + J{x)a ■ S 



ipQ{X) 



^ VQn{x)<pr^ 



[X) 



(12) 



n>l 



and 



[E - En - K-,- Un{x)] ipn{x) = Vnm{x)(Pm{x) ■ 



(13) 



Furthermore, relying on the large energy separation between the subbands in comparison 
with the magnitudes of Vnm{x) and J{x), we restrict our analysis to a two-subband model, 
keeping only the subband in the fixed wire whose energy is the closest to that of the swept 
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wire (the wave function of this subband is (fiix)) with the resulting set of equations given 

by 



E-Eo-K^- Uo{x) + JixYa ■ S 



[E-Ei-K^- Ui{x)] ^i{x) = Vix)Mx), 

where we have introduced V{x) = Voi(x) = Vio(x). 

Eqs. (14,15) can be decoupled using Green's functions: 



Ga(e) 



e — Kx — Uo{x) + J{x)a ■ S 



and 



Gi{e) = [e~Kx-U^{x)]-\ 
With these Green's functions Eqs. (14,15) can be formally integrated as 

0o{x) = G^{E - E^)V{x)iPi{x) 



and 



Accordingly, we obtain 



^^{x) = Gr{E ~ E{)V{x)^o{x). 



E-Eo-Kx~ Uo{x) + J{x)a ■ S 



and 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



^o(x) = V{x)Gi{E - Ei)V(x)^o(x), (20) 



[E-E^-Kx- f/i(x)] ^i(x) = V{x)Go{E - E,)V{x)^r{x). (21) 

The Green's function G'i(e) is a scalar Green's function, i.e. it is a unit matrix in the 
uncoupled spin space, whereas Go^e) has a more complicated structure. Nevertheless, it can 
be expressed in terms of two scalar Green's functions as 



Go(e) = -[3<7*(e)+/(e)l/ + 



a ■ S, 



where 



and 



g\e) = [e-Kx- U{x) + J{x)]-^ 
g'{e) = [e-K,-U{x)-?>J{x)]-' 



(22) 

(23) 
(24) 
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Now we are able to redefine the scalar potentials, as 



and 



f/o(x, E) = Uo{x) + V{x)Gi{E - E,)V{x) 



Uiix,E) = Uiix) + Vix)- [3g\E - Eo) + g%E - Eo)\ 



(25) 



(26) 



and introduce the tunneling-induced exchange coupling of electrons in the fixed wire to the 
local magnetic moment, 

1 r 



j{x, E) = -V{x)-^ [g\E - Eo) - g'{E - Eo)\ V{x). 



(27) 



As a result, we obtain equations for the description of electron dynamics in the swept and 
fixed wires in the form 



and 



E-Eo-K^- Uo{x) + J{x)a ■ S 



E-E^-K,-Ui{x)+j{x)a-S 



(po{x) 



ipi{x) 



0, 



0. 



(28) 



(29) 



Although the form of these two equations is very similar, and they can be both treated 
in the same manner ( as is discussed in the next Section), the results they yield will differ, 
depending on the specific shapes of the potentials and exchange couplings. In particular, 
while the shape of the coupling J{x) in Eq. PH|) is smooth, similar to that of the potential 
U{x), the exchange constant j{x) of Eq. is proportional to the potential V{x), and 

therefore is sharper than the bottleneck potential U{x). 



III. ELECTRON SCATTERING BY A LOCALIZED SPIN 

In this Section we determine the transmission coefficient and, correspondingly, the con- 
ductance of the quantum wire channel in the presence of electron scattering from a magnetic 
moment. The equation for the electron wave function in the uncoupled representation has 
the form 

[e-K^- U{x) + J{x)a ■ S) 0{x) = 0. (30) 

We are looking for a scattering solution of this equation with incident wave of the form 
Ae^''^, i.e. for a wave incident from x = — oo and having momentum hk = v2me. Here, A 
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is also the vector in the spin space, 



(31) 



with Wa = l^aP being the probabihty to have the certain initial orientation for the electron 
and magnetic moment spins. 

In the present paper we generalize the approach of Ref. (where the coupling constant 
J(x) was assumed to be a (5-function) to the case of a real spatial dependence of the coupling 
constant. To accomplish this, we find a solution of Eq. (j3(J|) for the wave function in the 
coupled representation and find the transmission coefficient for the wave function in the 
uncoupled representation by means of unitary transformation of Eq. (0). For the triplet 
and singlet states of the coupled representation, Eq. (jHUj) has the form 



-f/(x) + J(x))0f (x 







(32) 



and 



{e-K^- U{x) - 3J(x)) 0f (x) = 



(33) 



where indices t and s denote the triplet and singlet solutions, respectively. The scattering 
solutions of these equations have the following asymptotic behavior ^ 



lt,S+i 



T' Akx 



X] 



't,s< 



X 



+ Rt,s+e X 



and 



-ikx 



+ Rt s-e'''''' , X 



T' pikx 



► +00 

► — oo 

+ 00 

— oo 



(34) 



(35) 



and the transmission coefficient of the wave functions in the uncoupled representation can 
be obtained jlO] as 

^Tt 
Tt 








\{Tt + Ts) \{Tt-T, 



(36) 



y \{Tt-Ts) UTt + Ts) ) 
Now we can use the Landauer-Biittiker formula Il5| to determine the conductance of 
the quantum wire: 

(37) 



2e^ 

G = — ^ |To/3| 



a,/3 



where Wf^ gives the probabihty of initial spin configuration. With the transmission matrix 
given by Eq. (jHBj) . the conductance becomes 

^""^ Vtl' K + ^2) + ^ {\Ttf + \Tsf) {ws + W4)j (38) 



2 

4 



and, with account for the normahzation of probabihty, Y,a=i = 1, it takes the form 



2e^ 



(39) 



IV. CONDUCTANCE OF THE SWEPT WIRE 

In this section we show how our model for the local magnetic moment reproduces the 
correct conductance behavior when the swept wire is pinched off. As the basis for our 
analysis we use the expression for the conductance, Eq. (j39j) . jointly with our assumption 
that the potential Uo{x) and the exchange coupling J{x) are smooth functions in comparison 
to the electron wavelength, X = 27c/k. It should be noticed that the potential Uq{x) in Eq. 
fl28|) is not very smooth due to the additional contribution of a sharp potential V{x) (see 
Eq. (23)), however our analysis is still valid as long as this additional contribution is small 
in comparison to J{x). 

Functions Uq{x) and J{x) can be expanded into series near their maxima. Since the two 
functions are smooth we can assume that they take their maximum values at the same point 
X = 0. These expansions are given by 

, ^ x"^ d'^Uoix) , mufjx'^ , 

Uo{x) = f/o(0) + - ' U=o = Umax ^ (40) 

and 

T/ N T/«\ x'^d'^J(x), ^ muj']x'^ 

J(X) = J(0) + y^^U=0 = ^max j— ■ (41) 

Equations ()32|33|1 . defining the transmission coefficients Tt and T^, contain the effective 
potentials U^{x) = Uoix) — J{x) and f/_(x) = Uq{x) + 3J(x). In view of Eqs. (j4LH4ip . these 
potentials can be treated as inverse parabolic near their top, as 

TT ( \ TT T rnuj'ix'^ . 

^+l*^y ^ max J max ^ \ ) 

and 

TTlUj'^ X?' 

?7L(x) = Umax ~l" '^Jmax ^ ) (43) 
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where uj_ = yufj — ujj, uj+ = ^ujfj + ScUj. 

The transmission coefficient for the inverse parabohc barrier u{x) = is given 

by 



1 + e 



-27rr? 



-1/2 



(44) 



where rj = e/hu, and the energy, e, is measured from the top of the barrier. Thus, for our 
situation we obtain 



t 



^ Umax ~l~ Jmax 



and 



C Umax 3 Jmax 



(45) 



(46) 



The most important feature of these transmission coefficients is that the transmission prob- 
abihty, |t(?7)|^, is very close to the step function. 

Now we are able to calculate the conductance using Eq. (p^) . We assume the equivalence 
of all initial spin orientations, i.e. Wa = 1/4, and the conductance through the swept wire 
takes the form 



2e2 



t 



f Umax ~l~ Jmax 



1 

+ 4 



t 



C Umax 3 Jmax 



(47) 



The step-like structure of the transmission probability causes the conductance to repro- 
duce the step-like behavior of 0.7-anomaly. In case of ferromagnetic coupling between the 
electrons and local magnetic moment, Jmax > our model gives an additional conductance 
step at 0.75 x 2e'^/h, as 



2e2 



0, 



if ^ ^ Umax Jmax J 



0.75, if f/„ 



J m.ax *\ 6 *\ Ury 



(48) 



if 6 ^ Umax ~t~ 3 



It is interesting to point out that for antiferromagnetic coupling , Jmax < 0, we obtain 
the a conductance step at 0.25 x 2e^/h, which has been observed in experiments 13] and 
density-functional simulations ,4], as 



2e2 



, if 6 < Umax 3 1 1/ max \ ; 

0.25, if Umax 3 |t/j7jQia;| < £ < Umax ~\~ \Jmax \ i 
1 5 if 6 > Umax ~l~ Jmax ■ 



(49) 



10 



V. CONCLUSION 



In conclusion, we have examined electron dynamics in coupled quantum wires under con- 
ditions where a local magnetic moment is formed in one of the wires. In our theoretical 
model, the single-particle Hamiltonian has been employed with account of the specific ge- 
ometry defining the structure as well as electron scattering on the local magnetic moment. 
We have derived the coupled set of equations for the wave functions of electrons in both 
wires and have been able to decouple them obtaining equations describing electron dynam- 
ics in each of the two wires. While our analysis of electron processes in the fixed wire will 
be published elsewhere we determine the transmission coefficient and conductance of 
the swept wire (having the local magnetic moment) and show that our description repro- 
duces the experimentally observed features such as additional plateaus in the conductance 
at 0.75 ■ 2e^/h (for the ferromagnetic coupling between electrons and the local magnetic 
moment) and at 0.25 ■ 2e^ /h (in the case of the antiferromagnetic coupling). 
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FIG. 1: Electron micrographs of the critical region of the device, indicating the two measurement 
configurations 

FIG. 2: The figure on the left is a schematic of the experimental structure. The confining potential 
associated with this structure is modeled as a sum of potentials defining the two wires, W{y), 
(center) and defining the tunnel channel between them, V{x,y), (right). 
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Puller, et al., Figure 2 of 2. 
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